Manifestation of the shape and edge effects in spin-resolved transport through 

graphene quantum dots 
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We report on theoretical studies of transport through graphene quantum dots weakly coupled to 
external ferromagnetic leads. The calculations are performed by exact diagonalization of a tight- 
binding Hamiltonian with finite Coulomb correlations for graphene sheet and by using the real-time 
diagrammatic technique in the sequential and cotunneling regimes. The emphasis is put on the role 
of graphene flake shape and spontaneous edge magnetization in transport characteristics, such as the 
differential conductance, tunneling magnetoresistance (TMR) and the shot noise. It is shown that 
for certain shapes of the graphene dots a negative differential conductance and nontrivial behavior 
of the TMR effect can occur. 

PACS numbers: 73.63.Kv, 73.23.Hk, 73.22.Pr, 85.75.-d 
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I. INTRODUCTION 

Since its discovery^ graphene has been attracting an 
increasing attention due to its exceptional physical prop- 
erties and also possibilities of various promising practical 
applications^— For example, owing to a very long spin 
diffusion length observed in graphene^ one may expect 
that graphene will play an important role in future molec- 
ular spintronics. Moreover, with the advent of new pow- 
erful experimental techniques, it is possible now to engi- 
neer and fabricate graphene structures of various shapes 
and sizes, ranging from sheets of large area to extremely 
small graphene flakes. The latter can in particular ex- 
hibit single-electron charging effects and, thus, behave as 
typical quantum dots - similar to quantum dots based on 
two-dimensional electron gas£~— 

In such small graphene flakes, the role of edges is much 
increased in comparison to large graphene sheets. It 
is also well known on theoretical grounds that zigzag 
edges of graphene nanostructures have large densities of 
states, which in the presence of strong enough on-site 
Coulomb repulsions can result in the appearance of edge 
magnetism^ - — Indeed, it has been confirmed experi- 
mentally by scanning tunneling spectroscopy measure- 
ments that graphene nanoribbons and quantum dots re- 
veal highly enhanced densities of states (DOS) at the 
zigzag- type fragments of their edgesi 13 ' 14 It is worth not- 
ing, that the problem of graphene/graphite's edges has 
been recently under intensive studies, as the carbon- 
based nanostructures can potentially be used in modern 
nanoelectronics, including also spintronics J^r— Very re- 
cently, it has been demonstrated experimentally that the 
edge DOSs in graphene nanostructures are spin-split^ 
Following this line, in an attempt to gain additional in- 
sights into the edge states, we suggest another approach 
to the problem, namely a visualization of the effect of 
magnetic edges by the analysis of Coulomb blockade 
spectra for graphene dots of different geometries. To 
reach this objective we study the transport properties of 



graphene quantum dots coupled to ferromagnetic leads. 

As already mentioned above, in this paper we focus on 
the limit of rather small graphene flakes, and address the 
transport properties of graphene quantum dots weakly 
coupled to external ferromagnetic leads. The Coulomb 
blockade phenomena become then relevant. In particu- 
lar, we study the effects related with the shape and edges 
of the graphene flakes on various spin-resolved transport 
properties of the system, including differential conduc- 
tance, tunnel magnetoresistance (TMR) and shot noise 
(Fano factor). 

The question whether or not edge states can be probed 
by electronic transport methods is still a matter of in- 
tensive discussion. On the one hand, the edge mag- 
netism is critically suppressed in the case of contacts 
of good (or even moderate) transparency^ - — On the 
other hand, however, the edge magnetism appears in iso- 
lated graphene flakes and also survives when the flakes 
are weakly coupled to electrodes i 11 ' 23 ' 27 Additionally, 
the edge states may be localized and therefore (very) 
weakly conducting, which makes them hardly accessible 
by transport measurements. Our studies show that some 
information on the edge states can be extracted from 
transport measurements in the limit of weak coupling 
between the graphene flakes and electrodes, and when 
lateral dimensions of the flakes are not too large (com- 
parable to the localization length of the edge states). 
The first assumption makes the energy spectra of the 
flakes rather independent of the coupling to electrodes, 
while the second one makes the edge states accessible in 
transport measurements (albeit the corresponding con- 
ductance can be rather small). The Coulomb blockade 
spectra provide then a kind of unique shape-specific 'fin- 
gerprints', which also contain some information on the 
edge states. 

The paper is organized as follows. First, the model 
as well as the computational method based on real-time 
diagrammatic technique are briefly outlined in Sec. II. 
Then, the numerical results arc presented and thoroughly 
discussed in Sec. III. Summary and final conclusions are 
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FIG. 1: (Color online) Scheme of a graphene quantum 
dot coupled to external ferromagnetic leads. We consider 
graphene flakes of three different shapes: circular (a), rect- 
angular (b) and rhombic (c). A particular dot is coupled to 
external leads with coupling strengths described by and 
rj. for the left and right leads. It is assumed that the sys- 
tem can be in two magnetic configurations: the parallel and 
antiparallel ones, as sketched in the figure. 



in Sec. IV. 



II. THEORETICAL FRAMEWORK 

The considered system is displayed in Fig. [TJ It con- 
sists of a graphene flake that is weakly coupled to ex- 
ternal ferromagnetic leads. The coupling strengths are 
described by and Tf, for the left and right leads. We 
consider graphene flakes of three different hypothetical 
shapes: circular, rectangular and rhombic. It is assumed 
that the system can be in two magnetic configurations: 
either parallel or antiparallel one, see Fig. [T] 

In order to find the energy levels as well as magnetic 
moments of graphene quantum dots (GQDs) we have per- 
formed exact diagonalization of the following mean-field 
Hamiltonian 



H = - E i z ' ct > (*> j'I + \ E A v 



(1) 



Here, tij are the hopping integrals, \i,cr) stand for n- 
electron orbitals at site i with spin er, A, i<T = U (nj CT — 
rii-o-) describes the Stoner splitting, and n ia arc the 
respective occupation numbers. The latter have been 
computed sclf-consistcntly by summing up the squared 
eigenvectors corresponding to the eigenvalues not greater 
than the Fermi energy. The hopping integrals ti j are 
assumed to be nonzero only for nearest neighbors, and 
the nearest neighbor hopping parameter t is set to be 
equal t = 2.7 eV. In turn, the Coulomb on-site repul- 
sion is assumed to be U — 1.2i (see e.g. Ref. [Hj]). All 
the GQDs we consider are of comparable area (~9 nm 2 ) 
and consist of 350-400 carbon atoms. Here we focus on 
the shapes with a relatively small number of zigzag type 
edge atoms, and consequently few quasi-degenerate edge 
states in the vicinity of the Dirac point for U=0 (for 
triangular and hexagonal structures with purely zigzag 
edges see Refs. (27ll2a |). 



Exact diagonalization of the Hamiltonian ([I) yields the 
eigenvalues ei a that have been then used as an input to 
the mean-field Hamiltonian for the graphene quantum 
dot 



GQD 



Ec (XT \ 2 

~2~ ( N ~ n o) 



(2) 



Here, Ec is the phcnomcnological charging energy of the 
dot, eia is the energy of the dot's orbital discrete level 
I for spin a, ni„ = dj^d^ denotes the particle number 
operator for the level I, N = ^i-Tiia, and no is the 
number of electrons in an electrically neutral quantum 
dot. 

The leads are modeled by Hamiltonian of noninteract- 
ing quasi-particles 



H 



Leads 



E 2. 

X=L,R k,cr 



(3) 



where c^ k creates a spin-cr electron with wave vector k 
in lead a and e a u_ a is the corresponding energy In turn, 
tunneling processes between the dot and the leads are 
described by 
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(4) 



with v a i denoting the hopping matrix element between 
the dot level I and the lead a. The broadening of the 
GQD's levels can be described by = 27r j o^|u Q ;| 2 , 
where is the spin-dependent density of states in the 
lead a for spin subband a. In the case of ferromagnetic 



leads this can be then written as, T a 



+(-) 



(1 ± p a )T ai , 



with T a i = (r^ ; + r Qi )/2 and p a being the spin polar- 
ization of lead a. Here, (r~ ; ) corresponds to the 
coupling to the majority (minority) spin band. In the 
following we assume, r Q ; = r Q =T/2 and p a = p. Fur- 
thermore, in numerical calculations, for all three different 
shapes of the dots, we have also assumed the charging en- 
ergy Ec = 0.15 eV and the coupling strength V = 0.002 
cV. The Ec value has been found from a scaling law 
(against QD-size) formulated in Ref. (29|. Incidentally, 
this scaling leads also to acceptable estimations of charg- 
ing energies in Ref. @H[- 

In order to reliably determine the transport proper- 
ties of Coulomb-blockade graphene quantum dots weakly 
coupled to external leads, we employ the real-time dia- 
grammatic technique i^2r— This technique relies on sys- 
tematic perturbation expansion of the reduced density 
matrix and the operators of interest in the dot-lead cou- 
pling strength T. The calculation proceeds with the de- 
termination of respective self-energy matrices W, which 
enables the evaluation of the elements of the reduced den- 
sity matrix of GQD by using the following master-like 
equation^ 



(Wp s 



(5) 
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Here, p st is the vector containing stationary probabilities, 
W is the modified self-energy matrix W so as to include 
the normalization of probabilities and x = ^ a labels the 
states of the graphene quantum dot. Having found the 
occupation probabilities, the current flowing through the 
system can be calculated using the following equation^ 
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^Tr{WV}, 



(0) 



where W 1 is the modified self-energy matrix W to ac- 
count for the number of electrons transferred through 
the system. 

By performing the perturbation expansion of the re- 
spective self-energies one is then able to calculate the 
current order by order in tunneling processes. In this 
paper we have included the first and second-order self- 
energies. The first order of expansion corresponds to se- 
quential tunneling, which dominates transport outside 
the Coulomb blockade regime, whereas the second-order 
self-energies describe the cotunncling processes. Cotun- 
neling processes occur through virtual states of the sys- 
tem and are dominant in the Coulomb blockade regime.— 
Thus, to properly describe the transport properties of the 
system in the full range of bias and gate voltages, it is 
of vital importance to include both the first and second 
order terms of the perturbation expansion. 

With the aid of the real-time diagrammatic technique, 
we have determined the behavior of the current /, dif- 
ferential conductance G and tunnel magnctorcsistancc 
TMR of graphene quantum dots in both the linear and 
nonlinear response regimes. In addition, we have also 
calculated the bias and gate voltage dependence of the 
shot noise S and the corresponding Fano factor FM- The 
Fano factor is defined as, F — 5/(2e|/|), and describes 
deviation of the shot noise from the Poissonian value, 
relevant for uncorrelated tunneling events^ The main 
results and their discussion arc presented in the sequel. 



III. RESULTS AND DISCUSSION 

Below we present numerical results on transport 
through graphene quantum dots of similar sizes (in terms 
of the number of atoms), but different shapes. For sim- 
plicity, we assume the Fermi level of the leads to be equal 
to the Fermi level of a neutral graphene. We show that 
the corresponding energy spectra are strongly dependent 
on the GQD geometry, and that the emerging magnetic 
moments are essentially localized at zigzag-like segments 
of the dots' edges. This leads to different transport char- 
acteristics of particular GQDs, as shown and discussed 
in the following. More specifically, we show that the size 
of blockade regions (blockade diamonds in the bias-gate 
voltage dependence of the differential conductance) can 
be used to gain some information about the edge states. 
As it is well known, the size of the diamonds is deter- 
mined by the Coulomb charging energy Ec and the level 
spacing. Therefore, the sequence of diamonds depends 
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FIG. 2: (Color online) The atomic structure (a) and energy 
spectrum (b) of the rhombic graphene quantum dot. There 
are no magnetic solutions irrespective of the value of U, in- 
stead a quite pronounced energy gap opens. 



on whether a given level is spin degenerate or not. The 
absence of spin degeneracy, in turn, implies the presence 
of magnetic states. Since the edge states are either in the 
center of the energy gap or close to it, the corresponding 
diamonds can be easily identified. 



A. Rhombic graphene dots 

Let us begin our considerations with the case of 
graphene flake of rhombic geometry with exclusively 
armchair- type edges. The atomic structure of the 
graphene dot and the corresponding energy spectrum 
measured from the Fermi level, Ep — 0, is shown in 
Fig. [2l As one can see in this figure, the rhombic 
graphene flake has neither low-energy localized states nor 
magnetic moments at the edges. Moreover, the energy 
levels of the dot are independent of the Coulomb param- 
eter U, which is a consequence of the absence of magne- 
tized states. As one can readily see, such a structure has 
therefore a quite pronounced energy gap at the Fermi 
level. Thus, the armchair-edge rhombic geometry may 
be suitable for engineering graphene nanostructures use- 
ful for field effect transistor devices (with a pronounced 
ON/OFF current ratio). 

The bias and gate voltage dependence of the differen- 
tial conductance in the parallel and antiparallel config- 
urations is shown in Figs. and[3](b). The central 
white region corresponds to zero excess electrons in the 
dot. This region is relatively large due to the large en- 
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FIG. 3: (Color online) The bias V and gate V a voltage depen- 
dence of the differential conductance in the parallel (a) and 
antiparallel (b) configurations as well as the resulting TMR 
effect (c) calculated for rhombic- like graphene quantum dot. 
The parameters are: the charging energy Ec — 0.15 eV, the 
coupling strength T = 2 meV, the thermal energy fcsT = 20 
meV and the leads' spin polarization p = 0.5. 



ergy gap. With sweeping the gate voltage, one shifts 
the position of the graphene dot levels, changing thus 
the number of electrons in the dot. Due to particle- hole 
symmetry the spectrum is symmetric with respect to the 
central diamond. The second diamond (with increasing 
gate voltage starting from V g = 0) is much smaller, when 
compared to the middle one, as it is determined only by 
the charging energy Eq. The next diamond, apart from 
the charging energy, also includes the level spacing be- 
tween the first and second levels (for positive or negative 
energies). Since each level is spin degenerate, every sec- 
ond diamond (to the left or to the right of the central 
diamond) is determined only by the charging energy Ec 
and therefore all of them are relatively small. The other 
diamonds include additionally the level spacing and are 
therefore larger. Thus, the level spacings and the energy 
gap can, in principle, be determined from the size of the 
corresponding blockade diamonds. 

Except for particular structure of the Coulomb dia- 
monds related with the energy spectrum of the dot, which 
is independent of the magnetic configuration of the de- 
vice, one can see that the differential conductance in the 
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FIG. 4: (Color online) The Fano factor in the parallel (a) 
and antiparallel (b) magnetic configuration as a function of 
the bias and gate voltages for rhombic-like graphene quantum 
dot. The parameters are the same as in Fig. [3] 



parallel configuration Gp is larger than that in the an- 
tiparallel configuration Gap- This is related with the 
asymmetry in the couplings to the spin majority and 
spin minority bands of the ferromagnets. For the parallel 
configuration, the majority-majority channel is most con- 
ducting, while in the antiparallel configuration there are 
two weakly conducting majority-minority channels. As a 
consequence, Gp > Gap- The difference between the sys- 
tem transport properties in these two magnetic configu- 
rations is described by the tunnel magnctoresistance, de- 
fined as,? 7 TMR = (R AP -Rp)/R P , with R P and R AP de- 
noting the resistance in the parallel and antiparallel mag- 
netic configurations of the device. The bias and gate volt- 
age dependence of TMR is shown in Fig. (3{c) . First of all, 
one can see that the TMR is always positive, as in typical 
spin-value quantum dot devices^ The TMR is particu- 
larly enhanced in the central Coulomb blockade diamond 
[black area in Fig. [3jc)], where the sequential tunneling 
is exponentially suppressed while transport takes place 
via cotunncling events. In this transport regime the dot 
is empty and the current is driven by elastic cotunnel- 
ing processes ! 32 ' 36 The TMR is then exactly given by the 
Julliere valued TMR Ju11 = 2p 2 /(l - p 2 ), which for the 
assumed parameters yields 2/3. In other Coulomb di- 
amonds, the inelastic cotunncling processes become rel- 
evant and the TMR is generally smaller than Jullicre's 
value. Similar behavior can be observed in the sequential 
tunneling regime, where the TMR is much smaller than 
TMR Ju11 . This is related with the fact that the informa- 
tion about magnetic configuration of the device is now 
transferred by uncorrclatcd sequential tunneling events 
and the non-equilibrium spin accumulation which builds 
up in the graphene quantum dot. This mechanism is 
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less effective than direct spin-conserving cotunncling be- 
tween the left and right lead, thus the TMR is smaller 
than Julliere's value. 

The Fano factors in the parallel and antiparallcl mag- 
netic configurations of the device are shown in Fig. 3) 
In both cases, the shot noise in the blockade regions, 
where cotunneling processes dominate, is generally super- 
Poissonian, i.e. the corresponding Fano factors are larger 
than unity, F > 1. This is related with bunching of 
inelastic cotunncling events 3 -^ and has already been ob- 
served experimentally in transport through other quan- 
tum dot systems . 39 ' 40 Interestingly, in the central dia- 
mond, where elastic cotunneling mediates the current, 
the Fano factors in both magnetic configurations ap- 
proach unity. This is due to the fact that the elastic 
cotunneling processes are uncorrelated in time and the 
shot noise is Poissonian. In the regions, where trans- 
port is dominated by sequential tunneling processes, the 
corresponding shot noise is rather sub-Poissonian with 
the relevant Fano factors smaller than 1, F < 1. This 
suppression of the noise is due to Coulomb correlations 
between consecutive sequential tunneling processes. We 
note, that the calculations include also the thermal noise, 
which is dominant in the small bias voltage regime. Ac- 
cordingly, when V tends to zero, the current and shot 
noise both tend to zero, and the noise is dominated by 
the thermal noise (nonzero also at V = 0). Thus, the cor- 
responding Fano factor diverges when V — > 0, as marked 
with a horizontal stripe in Fig. |4j 



B. Circular graphene dots 
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FIG. 5: (Color online) The atomic structure (a) and energy 
spectrum (b) for circular graphene flake. In the case of U = 0, 
there are no magnetic edges, while magnetic edges appear for 
finite U. Note that in the latter case an energy gap (AE) 
opens. The corresponding magnetic configuration is also dis- 
played in (a). The maximum edge magnetic moments are 
equal to ca. 1/5 fiB- However, the total magnetic moment of 
the dot is equal to 1/j.b ■ 



Let us now have a closer look at the circular graphene 
flake. In this geometry the edge atoms have rather 
short zigzag-like (and armchair-like) coordinations. Con- 
sequently, any complete compensation of the edge mag- 
netization (antiferromagnetic alignment) can be hardly 
realized on geometrical grounds. Our calculations show 
that indeed the ground state magnetic configuration of 
the edges is ferromagnetic-like (fcrrimagnetic) , but with 
some admixture of antiparallel magnetic moments. In 
this situation there is no way to control spin dependent 
electric current with external magnetic field, unless one 
makes use of ferromagnetic electrodes. 

When U = (nonmagnetic dot), there is no gap at the 
Fermi level. However, as shown in Fig. [51 the presence 
of edge magnetism (U > 0) leads to the opening of an 
energy gap AE at the Fermi level. Noteworthy, near 
the half-band filling the highest occupied energy levels 
are fully spin-polarized - the dot behaves as a magnetic 
one. The overall magnetic moment of the dot shown in 
Fig. [5J is equal to 1/i_b. This, in turn, affects the spin- 
resolved transport properties of the system, which now 
show some asymmetry with respect to the bias and gate 
voltage reversal. 

Incidentally, in accordance with Lieb's theorem^ one 
can increase the net magnetic moment up to (Na — 



Nb)hb, by modifying the edges so as to increase the 
imbalance of the numbers Na and Nb of the A(B)- 
sublattice atoms. Our results are rather robust against 
moderate edge disorder. In particular, when a few edge 
atoms are removed then possible implications are consis- 
tent with predictions based on the aforementioned the- 
orem. Noteworthy, the asymmetry seen in Fig. [5] pro- 
vides direct visualization of the presence of GQD's un- 
compensated spin. In another context as the present 
one, the edge state effects on the electronic structure as 
well as charge and spin transport have been studied in 
Rcfs. gHH. 

The differential conductance spectra in both magnetic 
configurations are shown in Fig. [SJ First, we note that 
in equilibrium all the states of negative energy are occu- 
pied and the dot has one unpaired spin-down electron. 
The addition of a new electron costs then the charging 
energy plus half the level spacing AE, see Fig. [5j Note 
that the added electron is a spin- up one. This corre- 
sponds to the large central diamond in Fig. [5] The sub- 
sequent electron occupies spin-down level and the corre- 
sponding addition energy includes the Coulomb charging 
energy plus the spacing between the first and second lev- 
els of positive energy. Since the latter spacing is mach 
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FIG. 6: (Color online) The bias and gate voltage dependence 
of the differential conductance in the parallel (a) and antipar- 
allel (b) configurations as well as the resulting TMR effect (c) 
calculated for circular graphene quantum dot. The parame- 
ters are the same as in Fig.[3]except for thermal energy which 
is k B T = 10 meV. 



smaller than AE, the corresponding diamond is smaller 
than the central one, but it is still larger than the dia- 
mond determined by Ec only. It indicates that this state 
is magnetic. The next electron occupies again the spin- 
down level, which is almost degenerate with the previous 
level. The corresponding diamond is determined practi- 
cally by the charging energy Ec and is therefore smaller 
than the preceding one. The following two electrons oc- 
cupy the two spin-up states (also almost degenerate), so 
the corresponding two diamonds are determined by the 
charging energy plus the level spacing and charging en- 
ergy, respectively. Similar scenario holds for higher gate 
voltages, when higher energy levels become occupied by 
electrons, as well as to negative gate voltages. 

It is also worth noting that now the spectra are signif- 
icantly different from the corresponding ones for U = 0, 
see Fig. [5l Since for U = there is no gap at the Fermi 
level, the central diamond should be then determined by 
the charging energy only and therefore should be small. 
The second diamond should be relatively large while the 
three subsequent diamonds should be small and deter- 
mined by the charging energy. Thus, the conductance 
spectra can be used to distinguish the situation with 




FIG. 7: (Color online) The Fano factor in the parallel (a) 
and antiparallel (b) magnetic configuration as a function of 
the bias and gate voltages for circular graphene quantum dot. 
The parameters are the same as in Fig. [5] 



nonzero U from that with U = 0. 

Apart from the typical Coulomb blockade diamonds 
described above, the transport characteristics display an 
asymmetry with respect to the bias reversal, see Fig. [6] 
This is due to the spin splitting of the GQD's levels. 
In addition, in certain transport regimes one can find 
a negative differential conductance, which is present in 
both magnetic configurations. This effect is basically re- 
lated with the fact that for certain transport voltages 
electrons participating in transport are mainly spin-down 
(minority-spin) ones and the current is thus decreased. 

The difference between the currents flowing through 
the system in the two magnetic configurations of the de- 
vice results in the TMR effect shown in Fig. Because 
now the energy eigenstates of the graphene dot are no 
longer spin degenerate, one can observe a nontrivial be- 
havior of the TMR as a function of both the bias and 
gate voltages. First of all, it can be seen that in certain 
transport regimes the TMR can take values much larger 
than those given by Julliere's model. Furthermore, there 
are transport regimes where TMR changes sign and be- 
comes negative. To explain this behavior we will refer 
to approximate formulas for sequential tunneling and co- 
tunncling currents. Suppose the temperature is much 
larger than the coupling strength T, but still smaller 
than the charging energy Ec, Ec > ksT ^> T. Then, 
transport is mainly determined by thermally-activated 
sequential tunneling. The current is thus proportional to, 
/soq ~ J2a r L r fl/( r L + r fl)- If transport occurs through 
the spin-up level of the dot, the current in the parallel 
configuration is, /p 0q ~ (1 +p)T/2, while for the antipar- 
allel configuration one finds, /™p ~ (1 ~P 2 )r/2, which 
yields for the TMR, TMR SCC > = p/(l -p). If, in turn, the 
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FIG. 8: (Color online) The same as in Fig. [5] but for the 
rectangular graphene quantum dot. Now the maximum edge 
magnetic moments are equal to ca. 1/3 fiB- 



current is related with tunneling through spin-down lev- 
els, one gets 7p eq - (l-p)T/2 and I^p ~ (1 ~P 2 )r/2 for 
the parallel and antiparallel configurations, respectively, 
and the TMR is now negative, TMR soq = -p/(l +p). 
On the other hand, if the temperature is comparable to 
the coupling strength, k B T ss T E c , the current is 
mainly due to cotunneling processes. Its dependence on 
the spin polarization of the leads can be found if one con- 
siders elastic cotunneling regime, where I ~ J2a^l^R- 
Then, if cotunneling occurs through spin-up levels, the 
currents in the two configurations are proportional to, 
J P ~ (1 +p) 2 r 2 , 7 A p ~ (1 -P 2 )r 2 , yielding TMR = 
2p/(l — p), which for the assumed parameters (p = 0.5) 
gives, TMR = 2. Note that now TMR = 3 x TMR Ju11 , 
i.e. one obtains maximum TMR three times larger than 
in the case of rhombic-like graphene flake discussed in the 
previous section. However, if cotunneling occurs through 
spin-down levels, the current in the antiparallel configu- 
ration is the same as above, while the current in the paral- 
lel configuration becomes, 2p ~ (1 — p) 2 T 2 , which results 
in TMR = —2p/(l + p). For the assumed parameters 
one then finds, TMR = -2/3 = — TMR Ju11 , i.e. the sys- 
tem exhibits a large-magnitude negative TMR. Because 
the calculations are performed for Eq > k B T > T, one 
finds that the TMR in the Coulomb diamonds is much 
enhanced compared to Julliere's value, with 2p/ (1 — p) > 
TMR > p/(l — p), if transport is due to spin-up states 
of the graphene dot. On the other hand, if spin-down 
states are active in transport a negative TMR occurs 
with, -p/(l + p) < TMR < -2p/(l + p). These val- 
ues of the TMR can be clearly seen in Fig. [BJ:. 
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FIG. 9: (Color online) The bias and gate voltage dependence 
of the differential conductance in the parallel (a) and antipar- 
allel (b) configurations as well as the resulting TMR effect 
(c) calculated for rectangular graphene quantum dot. The 
parameters are the same as in Fig. [6] 

The corresponding Fano factors in both magnetic con- 
figurations are shown in Fig. [7J In both cases, the shot 
noise is generally super- Poissonian in the Coulomb block- 
ade regions, where cotunneling processes dominate, and 
bunching of inelastic events enhances the noise. In other 
regimes transport is dominated by sequential tunnel- 
ing processes and the corresponding shot noise is sub- 
Poissonian, F < 1. As before, this is due to the Coulomb 
correlations between sequential tunneling events and the 
Pauli exclusion principle, resulting in suppressing of the 
noise. 



C. Rectangular graphene dots 

In this subsection, transport through a rectangular 
graphene quantum dot with zigzag and armchair edges is 
studied. In contrast to the cases discussed earlier, with 
cither nonmagnetic edges (rhombus dot) or nearly fer- 
romagnetically aligned ones (circular dot), the present 
case corresponds to a spontaneous antiparallel alignment 
of magnetic edges. This device may therefore act as an 
effective spin- valve or a magnetoresistive sensor, and ap- 
pears potentially useful for spintronics i 16 ' 34 
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FIG. 10: (Color online) The Fano factor in the parallel (a) 
and antiparallel (b) magnetic configuration as a function of 
the bias and gate voltages for rectangular graphene quantum 
dot. The parameters are the same as in Fig. [S] 

Figure |S] presents atomic structure of the considered 
graphene flake and the corresponding energy spectrum. 
It can bee seen that for a sufficiently large on-site 
Coulomb repulsion parameter U, an energy gap opens in 
the energy spectrum, resulting in a metal/semiconductor 
(zero-gap/finitc-gap semiconductor) transition. More- 
over, the energy levels of this graphene dot remain spin 
degenerate (perfect compensation of the net magnetiza- 
tion), in contrast to the case of the circular dot. 

The subsequent figures show the corresponding bias 
and gate voltage dependence of the differential conduc- 
tance and TMR, see Fig. [9l as well as the Fano factor, 
sec Fig. [TUJ Since there is a gap at the Fermi level, the 
central diamond in the differential conductance is now 
rather large and includes the charging energy plus a half 
of the energy gap. Because the energy levels are now 
spin degenerate, every second diamond is smaller, as it 
is determined merely by the charging energy Eq- More- 
over, because of the spin degeneracy of the levels, the 
differential conductance in the parallel configuration is 
always larger than that in the antiparallel configuration, 
yielding positive TMR, sec Fig. M The behavior of TMR 
resembles that in the case of the rhombic graphene flake. 
The TMR is given by Julliere's value in the Coulomb 
diamonds where elastic cotunneling drives the current, 
and is much more suppressed in the other diamonds with 
either inelastic cotunneling or sequential tunneling. 

The bias and gate voltage dependence of the Fano 
factor in both magnetic configurations of the device is 
shown in Fig. [TO] The general features of these spec- 
tra are qualitatively similar to those discussed earlier in 
the case of the rhombic-like graphene dots: F « 1 in 
the Coulomb blockade regions when elastic cotunneling 



is relevant, F > 1 in the blockade regions when inelas- 
tic cotunneling is present, and F < 1 in the sequential 
tunneling regime. 



IV. CONCLUSIONS 

In this paper we have analyzed transport properties 
of graphene quantum dots of different geometries in the 
presence of the on-site Coulomb repulsion. Characteris- 
tic features due to GQD edges are shown to be reflected 
in the Coulomb blockade diamonds and other transport 
characteristics. It has been demonstrated that the TMR 
effect in graphene quantum dots may be quite large, 
much larger then the corresponding Julliere's value of 
TMR. In turn, the shot noise was found to take both sub- 
Poissonian as well as super- Poissonian values. Out of the 
three graphene flake shapes studied here, the armchair- 
edge rhombus flake shows no magnetism, the circular 
flake has uncompensated ferrimagnetic edge configura- 
tion, whereas the rectangular flake displays antiparallelly 
aligned zigzag-edge magnetizations. All these features 
manifest themselves in transport characteristics, provid- 
ing information about the geometry and edge states of 
graphene dots. 

In frame of the model assumed, the conductance spec- 
tra can be used to check whether the energy spectrum 
includes zero energy levels or not. If they do, then the 
first diamonds are determined only by the charging en- 
ergy and are small. If not, then the central diamonds are 
large as they include not only the charging energy, but 
also the level separation due to the size quantization. 
However, the model is not free from some limitations, so 
the conductance spectra in real situations may be differ- 
ent. Accordingly, the conclusions concerning magnetic 
edge states should rather be regarded as indicative only. 
First of all, the description is based on mean-field approx- 
imation for charging energy Eq and also for the coupling 
parameter T. In reality, both parameters may depend on 
the particular states taking part in transport. As already 
mentioned in the introduction, this may be especially im- 
portant for r, when transport occurs via localized edge 
states. As the variation of Ec has some influence on the 
size of the Coulomb diamonds, the variation of T mod- 
ifies mainly the conductance and TMR, leaving the size 
of the diamonds unchanged. Being aware of all the un- 
derlying assumptions, we believe our results present a 
sound starting point aimed at transport characterization 
of the electronic states (including also the edge states) of 
graphene flakes. Since both the crucial parameters (es- 
pecially the coupling parameter T) depend on a particu- 
lar experimental realization, a more accurate description 
can be performed when the relevant experimental data 
are available. 
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